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The ability to break symmetry and polarize through 
self-organization is a fundamental feature of cellu- 
lar systems. A prevailing theory in yeast posits that 
symmetry breaking occurs via a positive feedback loop, 
wherein the adaptor protein Beml promotes local activa- 
tion and accumulation of Cdc42 by directly tethering 
Cdc42'^^'' with its guanine nucleotide exchange factor (GEF) 
Cdc24. In this paper, we find that neither Beml nor the 
ability of Beml to bind Cdc42*^^^ is required for cell polar- 
ization. Instead, Beml functions primarily by boosting GEF 



activity, a role critical for polarization without actin fila- 
ments. In the absence of actin-based transport, polariza- 
tion of Cdc42 is accomplished through Rdil, the Cdc42 
guanine nucleotide dissociation inhibitor. A mathematical 
model is constructed describing cell polarization as a 
product of distinct pathways controlling Cdc42 activation 
and protein localization. The model predicts a nonmono- 
tonic dependence of cell polarization on the concentration 
of Rdil relative to that of Cdc42. 



Introduction 

Actin and microtubule cytoskeletons are key players driving 
cellular symmetry breaking in diverse morphogenetic systems 
(Li and Gundersen, 2008). Cytoskeletal polymers, through their 
inherent structural and biochemical polarity, orchestrate movement 
and localization of regulatory molecules, which in turn govern 
cytoskeleton assembly and organization. Such positive feedback 
loops between structural and regulatory components are features 
in most models of cell polarization (Wedlich-Soldner et al., 2003; 
Marco et al., 2007; Seetapun and Odde, 2010; Goehring et al., 
2011; Ku et al., 2012; Jose et al, 2013; Slaughter et al., 2013). 
The budding yeast, however, presents a rather unique case in 
which cellular symmetry breaking can be accomplished with or 
without contribution from the cytoskeleton (Ayscough et al., 
1997; Irazoqui et al., 2003; Wedlich-Soldner et al., 2004). 

Yeast cells naturally bud in a defined pattern, termed bud 
site selection, guided by sites of previous cell division — bud 
scars (Casamayor and Snyder, 2002). However, cells are also 
known to polarize with equal efficiency, albeit in random orien- 
tations, when bud scar cues are eliminated or cells lose the abil- 
ity to recognize them, reflecting a self-organizational mechanism 
at work. A highly conserved regulator of cellular symmetry 
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breaking and polarized morphogenesis is the small GTPase 
Cdc42 (Etienne-Manneville, 2004). Localization of Cdc42 from 
an isotropic distribution around the cortex to a focused cap on 
one side of the cell is thought to be the key step in symmetry 
breaking, orchestrating, through Cdc42 effectors, a massive 
polarized reorganization of cellular components such as the actin 
cytoskeleton and the secretory machinery, leading to the initia- 
tion of polarized growth to generate the bud. The Cdc42 polar 
cap is dynamically maintained: diffusion of Cdc42 within the 
membrane and exchange with the fast-diffusing cytosolic pool 
must be balanced by continuous recycling and retargeting of 
Cdc42 to the cap region (Wedlich-Soldner et al., 2004; Marco 
et al., 2007; Slaughter et al., 2009). 

Two distinct models have been proposed to explain sym- 
metry breaking in yeast. One model relies on a positive feedback 
loop between Cdc42-regulated assembly of the actin cytoskele- 
ton and actin cable-mediated transport of exocytic vesicles car- 
rying Cdc42 to the nascent cap (Wedlich-Soldner et al., 2003; 
Marco et al., 2007). Crucial to this model, recent studies found 
that membrane microdomains and endocytic corralling serve to 
maintain Cdc42 concentration in the polar cap (Jose et al., 2013; 
Slaughter et al., 2013). In normally growing cell populations, 
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Figure 1 . Bemi and GEF localization ore not required for symmetry breaking. (A) Cells with indicated genotypes were plated on YPD media and grown 
at 23°C for 3 d. Cells were grown overnight in liquid culture and then diluted to an OD of 1 . This and a series of 1 0-fold dilution were spotted left to right. 
(B) Polarization of GFP-Cdc42 wild type, axl2A raxTA, and rsrlA cells after release from Gl arrest into media containing DMSO or LatA. The percentage 
of cells with polarized GFP-Cdc42 was determined at different time points (given in minutes] after release. The plots show means from averaging three 
experiments, and error bars correspond to SEMs. More than 80 cells were counted per time point per experiment. (C) Maximum projections of represen- 
tative polarized cells from A. (D) Representative kymographs of polarizing GFP-Cdc42 in axl2A raxlA and rsrlA cells in LatA. Note the unstable polar 
cap in axl2A raxlA relative to that in rsrlA. See also Fig. SI and Videos 1 and 2. (E) Quantification of cap duration in LatA of polarized GFP-Cdc42 in 



1 DBS JCB • VOLUME SOS • NUMBER 7 • 201 3 



however, disruption of actin only renders polarization less 
efficient (Wedlich-Soldner et al., 2004). An actin-independent 
model for cell polarization is centered on Cdc24, the lone Cdc42 
guanine nucleotide exchange factor (GEF) in yeast required for 
converting Cdc42 into the active, GTP-bound form. Cdc24 has 
been proposed to form a complex with Beml (Peterson et al., 
1994; Zheng et al., 1995; Ito et al., 2001), an adaptor-like protein 
sharing several protein interaction motifs (FBI, SH3, and PX) 
with animal cell polarity protein PAR6 and the p67Phox adaptor 
protein in the NADPH oxidase complex. Beml also has the abil- 
ity to bind Cdc42'''^'' (Bose et al., 2001; Yamaguchi et al., 2007), 
like PAR6 (McCaffrey and Macara, 2009), as well as effectors 
of Cdc42 such as the p21-activated kinase Cla4 (Bose et al., 
2001). The model posits that positive feedback occurs when the 
adaptor-GEF complex is recruited to an initial small accumula- 
tion of Cdc42'^^'' through the interaction with Beml, leading to 
GEF localization and thus autocatalytic Cdc42 activation at the 
nascent cap (Butty et al., 2002; Kozubowski et al., 2008). Beml 
was found to be essential for polarization only in the presence of 
Latrunculin A (LatA), leading to the proposal that Beml and 
actin represent two parallel pathways (Wedlich-Soldner et al., 

2004) . However, several more recent studies proposed the 
Beml -mediated positive feedback loop to be the sole mecha- 
nism for symmetry breaking in yeast on the basis of synthetic 
lethality between BEMl and RSRl, encoding a Ras-like GTPase 
required for bud site selection (Kozubowski et al., 2008; Howell 
et al., 2009). This synthetic lethality was interpreted as to show 
that Beml is required for cell polarization in the absence of bud 
site selection. 

Another unresolved question is how an autocatalytic acti- 
vation of Cdc42 might lead to polarized localization of Cdc42. 
Soluble Cdc42 in the cytosol is chaperoned by Rdil (Koch 
et al., 1997; Slaughter et al., 2009; Das et al., 2012). Rdil, like 
other Rho GDP dissociation inhibitor proteins, has the ability to 
extract Cdc42 from the cortex (DerMardirossian and Bokoch, 

2005) , and this activity is facilitated by modulation of electro- 
static interactions of Cdc42 with lipids by the lipid flippase 
Lem3 (Das et al., 2012). The Rdil-mediated Cdc42 recycling is 
required for maintaining the polar cap when actin is inhibited 
with LatA (Slaughter et al., 2009; Das et al., 2012). Rdil may 
passively facilitate polarization if it is prevented from extract- 
ing Cdc42 at the cap through a GEF-dependent mechanism 
(Kozubowski et al., 2008; Savage et al., 2012), but the rapid. 



Rdil -dependent exchange of Cdc42 at the polar cap (Wedlich- 
Soldner et al., 2004; Slaughter et al., 2009) argues against this 
idea. In this work, we critically examine the role of Beml in 
symmetry breaking with or without actin and present an alterna- 
tive model for Cdc42 polarization centered on Rdil-mediated 
Cdc42 recycling. 

Results 

Bem1 is not required for cell polarization 
in tl-ie absence of tlte bud scar cue 

Axial and bipolar budding patterns rely on the trans-membrane 
markers Axl2 and Raxl-Rax2 complex, respectively (Roemer 
et al., 1996; Sanders et al., 1999; Kozminski et al., 2003; Fujita 
et al., 2004; Gao et al., 2007). Transmission of the bud scar sig- 
nal to Cdc42 requires Rsrl (Kozminski et al., 2003; Park and 
Bi, 2007), forming the basis of using rsrl A as the genetic back- 
ground for studying spontaneous symmetry breaking. However, 
Rsrl may have a more general role in polarization: Rsrl physi- 
cally interacts with both Cdc42 (Kozminski et al., 2003) and 
Cdc24 (Zheng et al., 1995; Park et al., 1997), as well as Beml 
(Park et al., 1997), and deletion of RSRl was reported to result 
in spatially unstable Cdc42 caps (Ozbudak et al., 2005; Howell 
et al., 2012), particularly in the presence of LatA. Thus, rsrlA 
could compromise both bud site selection and the core mecha- 
nism of cellular symmetry breaking. 

To create a more specific genetic background for studying 
cellular symmetry breaking, we deleted the bud scar markers 
Axl2 and Raxl to remove the structural cues for both axial and 
bipolar budding, as this leaves the Rsrl GTPase module intact. 
We confirmed that axl2A raxlA cells bud in random orienta- 
tions (Fig. SI A). Like rsrl A cells, axl2A raxl A cells grew 
similarly to the wild type under the standard experimental con- 
dition (Fig. 1 A). We next measured the rate of polarization on 
the population level in the presence or absence of LatA after 
release from a pheromone-induced Gl arrest. axl2A raxl A cells 
polarized at rates similar to wild-type cells in both DMSO and 
LatA. In contrast, rsrl A cells showed slower polarization in the 
presence of LatA (Fig. 1, B and C; and Fig. SI B). Time-lapse 
imaging was then performed to examine polarization of GFP- 
Cdc42 in single polarizing rsrl A or axl2A raxl A cells in LatA. 
The single-cell analysis revealed a notable difference in polar 
cap stability between the different mutants: polar caps in rsrlA 



axl2A raxl A and rsrl A cells from videos such as the ones shown in D and E. The maximum intensity of the polar cops was measured and plotted over 
time and was fitted to Gaussian curves. The duration of the cap is reported as the full width half maximum of the Gaussian fit. Boxes, SEM; whiskers, SD; 
horizontal lines, median values. P < 10"*. For each dataset, n > 20 cells taken from three experiments. (F and G) Time-lapse imaging of Cdc24-GFP in 
mid-log phase axl2A raxlA bemlA (F) and rsrlA bemlA (G) cells. Notice the lock of Cdc24 polarization in G despite successful cell polarization and 
budding compared with strong Cdc24 polarization in F. Arrows point to incipient bud sites. White outlines indicate the perimeter of the budded cell. See 
also Videos 3 and 4. (H] Kymographs of polarizing cells from F and G. Shown at the bottom ore intensity profiles token from the time point immediately 
before bud emergence, indicated in the kymograph by arrows. (I) Quantification of polarization of Cdc24-GFP in axl2A raxl A bem lA and rsrl A bem lA 
cells before bud emergence. Intensity at the emergent bud site is normalized to cortical intensity of the rear half of the cell. Three data points are shown 
for each cell, corresponding to time points 2, 4, and 6 min before bud emergence, n = 1 1 cells for axl2A raxlA berr)lA and 17 cells for rsrlA berr)lA. 
Boxes, SEM; whiskers, SD; horizontal lines, median values. P < 10"^. (J) Maximum projections of representative cells with polarized Cdc24-GFP in ax/24 
raxl A and rsrl A cells in LatA or DMSO. (K) Quantification of polarization of Cdc24-GFP in axl2A raxl A and rsrl A cells in LatA 50 min after pheromone 
arrest release from maximum projections of z-stack images. Cells with polarized GFP-Cdc24 were identified, and cortex masks were generated from a 
separate fluorescent channel (mCherry-Cdc42; see Materials and methods]. Mean intensity of the cortex within the mask was determined in 1 0° increments 
centered on the polar cap and normalized to the rear half as in I. Normalized profiles were averaged over n > 1 8 cells. Error bars correspond to SEMs. 
For comparison of values at 0°, P < 0.00 1 . Bars, 5 \im. 
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Figure 2. Bern 1 -Cdc42*'"' binding Is not required for GEF localization or cell polarization. (A) Schematic drawing of the proposed Berrl feedbacl< loop. 
(B) Beml domains and l<nown interacting partners. (C) Bern! mutations used in this study and their projected effects on Beml interactions. (D) Polarization 
of GFP-Cdc42 in axl2A rax ]A BEMI , axl2A rax lA bem lA, and ox/2zl rax lA bem cells, as in Fig. 1 B. The plots show means and SEMs from three 
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lasted for 440 ± 60 s (n = 24) and then sometimes reappeared at 
a different location (Fig. 1, D and E; Fig. SI C; and Video 1), 
whereas axl2A raxlA cells maintained stable polar caps for 
much longer periods, averaging 1,720 ± 170 s (n = 20; Fig. 1, D 
and E; Fig. SID; and Video 2). This result shows that Rsrl has 
a function in symmetry breaking in addition to its role in bud 
site selection. 

Although extensive data were interpreted based on the 
assumption that Beml is required for viability in the rsrl A 
background (Kozubowski et al., 2008; Howell et al., 2009, 
2012), both beml A rsrl A double and beml A axl2A raxlA tri- 
ple mutants are viable and grew similarly to bemlA alone in 
the S288c genetic background (Fig. 1 A and Fig. SI E), sug- 
gesting that Beml is nonessential for symmetry breaking. 
Strikingly, time-lapse imaging revealed that bemlA rsrlA cells 
were able to polarize and bud with severely reduced or unob- 
servable Cdc24 localization to the incipient bud site (Fig. 1, 
G-I; and Video 3). In contrast, beml A axl2A raxlA cells showed 
prominent polarization of Cdc24-GFP before bud emergence 
(Fig. 1, F, H, and I; and Video 4), similar to beml A cells as 
previously reported (Gulli et al., 2000; Butty et al., 2002). 
Polar localization of Cdc24-GFP in LatA was also significantly 
reduced in rsrlA cells compared with axl2A raxlA cells 
(Fig. 1, J and K). Collectively, these results show that, although 
Rsrl and Beml together control Cdc24 localization in the polar 
cap, cells could still polarize and bud without local concentra- 
tion of the GEF. 

Symmetry breaking does not rely on Beml- 
mediated Cdc4S-t:a-CdcS4 feedback loop 

All existing models assume that Beml functions in symmetry 
breaking by mediating a positive feedback loop connecting 
Cdc42 and Cdc24 as outlined in Fig. 2 A. Whereas this assump- 
tion was supported by gain-of-function experiments in which, 
for example, Cdc24 covalently linked to Cla4 could bypass 
the requirement for Beml for viability in rsrl A background 
(Kozubowski et al., 2008), the assumption has not been rigor- 
ously tested using loss-of-function approaches under more 
physiological settings. To this end, we pursued an unbiased 
investigation into how Beml participates in symmetry breaking 
by systematically disrupting each of the known physical inter- 
actions of Beml (Fig. 2 B) using specific point mutations vali- 
dated in previous studies (Fig. 2 C; Butty et al., 2002; van 
Drogen-Petit et al., 2004; Yamaguchi et al., 2007). We note that 
all beml mutants analyzed in the following experiments were 
in the axl2A raxlA background for observing symmetry break- 
ing in the absence of the bud scar cue. 

We began by disrupting the binding of Beml to Cdc42'^'^^, 
a critical step in the proposed Beml -dependent feedback loop. 



This was achieved by replacing BEMl with a centromeric plas- 
mid expressing beml (under the BEMl promoter) bearing the 
N253D mutation, which lies in the Cdc42 interaction domain 
and was previously shown to abolish Beml binding to active 
Cdc42 (Yamaguchi et al., 2007). Surprisingly, polarization of 
GFP-Cdc42 was unaffected in beml'""^'^ cells in the axl2A 
raxlA bemlA background compared with BEMl cells in the 
same background in the absence or presence of LatA (Fig. 2, 
D and E). Cdc24-GFP localization was reduced in beml'^'^^" 
cells in the axl2A raxlA bemlA background in LatA (Fig. 2, 
F and G) but remained higher than in rsrl A cells (compare Fig. 2 G 
with Fig. 1 K). If the N253D mutation indeed disrupted the 
Beml-Cdc42 interaction, Beml^^^^^'-GEP should be unable to 
polarize to the Cdc42 cap. Indeed Beml'^^'^'-GFP polarized 
poorly (Fig. 2 H), and we quantified the polarization strength of 
BemP'^^^-GFP compared with Beml-GFP, which indicated a 
significant (P < 10"*'*) reduction of polarization of Beml''*^^'°- 
GFP compared with Beml-GFP (Fig. 2 I). This supports the 
in vivo disruption of the Cdc42-Beml interaction and shows 
that symmetry breaking does not require localization of Beml 
to the site of Cdc42 accumulation. 

Bem1-CdcS4 interaction contributes to 
polarization by boosting GEF activity 

We next disrupted the interaction between Beml and Cdc24 
mediated through heterodimerization between PBl domains 
located at the C termini of these proteins (Ito et al., 2001; Butty 
et al., 2002). Supporting a crucial role for this interaction in 
symmetry breaking, mutations in Beml abolishing PBl domain 
binding (beml'^^'^"^'^^'''^'') resulted in failure of cells to polarize in 
the presence of LatA, whereas polarization in these mutants 
was only reduced without LatA, similar to that in bemlA 
(Fig. 3, A and B). 

Because Cdc24 was able to localize even in the absence of 
Beml localization (Fig. 2 F), it was unlikely that the failed po- 
larization of PB 1 domain mutant cells was caused by a lack of 
Cdc24 localization. A previous study hypothesized that Cdc24 
exists in an autoinhibited form and that binding by Beml might 
help to relieve this inhibition (Shimada et al., 2004). To investi- 
gate this possibility, we developed a Forster resonance energy 
transfer (FRET)-based biosensor (Itoh et al., 2002; Hodgson 
et al., 2008) for Cdc42 activation that consists of a linked con- 
struct of yeast Cdc42 with the Cdc42/Rac-interactive binding 
(CRIB) domain from Cla4 (Fig. 4 A), which interacts only with 
the active, GTP-bound form of Cdc42. These are flanked by 
GFP and mCherry, such that when the CRIB and Cdc42°'^^ 
within the sensor are bound, the GFP and mCherry are brought 
in close proximity for energy transfer to occur, as illustrated in 
Fig. 4 A. The polybasic-CAAX box region of Cdc42 was moved 



experiments, with n > 80 cells counted per time point per experiment. (E) Maximum projections of representative polarized cells from D. (F) Maximum 
projections of representative cells with polarized Cdc24-GFP in axl2A rax TA bem Z'^^^^" cells after pheromone arrest release in LatA or DMSO. (G) Polar- 
ization of Cdc24-GFP in axl2A rax lA cells and in axl2A rax lA bem /'^^■'^^ cells in LatA, as in Fig. 1 K. For the comparison of peak values at 0°, P < 0.01 . 
(H) Localization of Beml-GFP or bem /'^^■'■"'-GFP in cells with polarized mCherry-Cdc42 after pheromone arrest release in LatA. Z-stack images of represen- 
tative cells were subjected to a 1 x 1 Gaussian filter before maximum projection. (I) Histogram of Beml-GFP polarity (maximum pixel intensity/mean pixel 
intensity per cell) for Beml-GFP and bem /'^^^""'-GFP at 50 min after pheromone arrest release in LatA. n = 520 cells were quantified for each genotype. 
Bars, 5 pm. 
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Figure 3. Beitil contributes to symmetry breaking by boosting the GEF activity of Cdc24. (A) Polarization of GFP-Cdc42 in axl2A rax/4 BEMl cells, 
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Experimental procedure was as described in Fig. 1 B. Tfie plots show means from averaging three experiments, and error bars correspond to SEMs. AAore 
than 80 cells were counted per time point per experiment. (B] Maximum projections of representative cells from A. (C) Localization of Cdc24'^™'-GFP 
in faem /'^''^''^'^^'''''' cells with polarized mCherry-Cdc42 after pheromone arrest release in LatA. Z-stack images of representative cells were subjected to a 



1 X 1 Gaussian filter before maximum projection. (D) Polarization of Cdc24'^™'-GFP and mCherry-Cdc42 in bem 
Fig. 1 K for n = 15 cells. Error bars show SEA/\s. Bars, 5 \im. 
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C terminally to mCherry to allow proper prenylation and mem- 
brane anchorage. Higher FRET efficiency indicates higher net 
GEF activity toward Cdc42. FRET was measured using the ac- 
ceptor photobleaching approach as described in the Materials 
and methods (Fig. 4, A-C; and Fig. S2). Mutations were intro- 
duced into the Cdc42 portion of the biosensor, resulting in 
either a constitutively GTP-bound (Q61L) or GDP-bound (D57Y) 
state to serve as positive or negative controls, respectively. The 
positive and the negative controls showed well- separated FRET 
efficiencies (n > 25 cells for each strain; P < 10"*), whereas the 
wild-type biosensor showed a FRET level intermediate between 
the controls as expected (Fig. 4 D). 

We introduced the aforementioned Cdc42 activation bio- 
sensor into various mutant strains. Deletion of RSRl resulted in 



FRET measurements similar to that in wild-type cells (P = 0.3), 
whereas deletion of BEMl resulted in a significant lower FRET 
efficiency (P < 10"'; Fig. 4 D). Interestingly, disruption of the 
Beml-Cdc24 interaction through the heml'^'*^'^'^^*'^'' mutation 
reduced the biosensor FRET by an extent similar to that by 
bemlA (P = 0.8 compared with bemlA, and P < 0.01 compared 
with BEMl; Fig. 4 D). In contrast, GEF activity in beml^"^° 
cells was not significantly different from the wild type (P = 0.4; 
Fig. 4 D). These results show that the Beml interaction with 
Cdc24 is primarily responsible for enhancement of the latter' s 
GEF activity. Because it was proposed that the PB 1 domain of 
Cdc24 acts as an autoinhibitory domain and that binding by Beml 
may help relieve the autoinhibition (Shimada et al., 2004), dele- 
tion of the PBl domain should result in a constitutively active 
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Figure 4. Biosensor measurements of Cdc24 GEF activity level in various strains. (A) Schematic representation of the FRET-based biosensor for Cdc42 
activation. Active Cdc42°^'' binds the CRIB domain, bringing the flanking GFP and mCherry into proximity and allowing FRET. (B) A representative wild- 
type cell expressing the biosensor containing wild-type Cdc42. Leftmost image shows the sum of the time series for GFP. The center image shows the same 
image with cortex mask applied (see Materials and methods). The right image shows the FRET efficiency at each pixel within the cell as indicated by the 
heat bar. Bar, 5 \im. (C) Mean FRET efficiency (orange curve) was measured and plotted in 10° increments along the cortex in the masked image (see B, 
center image). Mean normalized GFP intensity at the cortex (blue curve) was also measured and plotted. Each plot shows means from >25 cells, and error 
bars show SEMs. (D) FRET efficiency in the cap region in indicated strains. The left two bars were from wild-type cells expressing positive and negative 
control sensors. Mean FRET efficiency was measured within the yellow region highlighted in C (see Materials and methods for details). Each histogram 
shows means from >25 cells, and error bars show SEMs. *, P < 0.01 ; **, P < 0.001 ; * **, P < 1 0"^. 



version of Cdc24. Indeed, Cdc24^''^' expressed under the 
CDC24 promoter from a centromeric plasmid boosted biosen- 
sor FRET in CDC24 beml^-^^"'^-''^''^'' cells (P = 0.1 compared 
widi BEMl , and P < 0.01 compared with beml^'*^'"^^^''^''; Fig. 4 D). 
Importantly, the deregulated GEF (CDC24^''^') partially res- 
cued cell polarization in beml^'^^''^'''^''^'' cells in the presence of 
LatA (Fig. 3, A and B), supporting the notion that Beml's role 
in cell polarization without actin is mainly mediated through its 
interaction with Cdc24, which stimulates Cdc24 GEF activity. 



Cdc24'^™'-GFP localized poorly to caps of mCherry-Cdc42 in 
LatA-treated cells (Fig. 3, C and D), further suggesting that cells 
can polarize without GEF localization if the overall cortical 
GEF level is high. 

A second mechanism by which Beml has been suggested 
to modulate Cdc24 activity is mediation of complex formation 
between Cdc24 and the Cdc42 effector Cla4. Cdc24 is hyper- 
phosphorylated by Cla4 in this complex, which requires direct 
binding of both Cla4 and Cdc24 by Beml (GuUi et al., 2000; 
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Bose et al., 2001; Wai et al., 2009). Expression of a fusion con- 
struct of Cdc24'^''^' linked with Cla4 was shown to rescue the 
synthetic lethality of Arsrl Abeml in the yeast strain back- 
ground used in a previous study (Kozubowski et al., 2008), 
which was interpreted as supporting a sufficiency of complex 
formation between Cdc24 and Cla4 for symmetry breaking. 
However, we found that the rescued cells failed to polarize 
in the presence of LatA (Fig. 5, A and B), suggesting that 
the fusion construct relies on an actin-dependent mechanism to 
accomplish the rescue. Because Beml and Cla4 interact through 
the Beml SH3 domain and a proline-rich motif in Cla4 (Bender 
et al., 1996), we tested the disruption of the projected Beml- 
Cla4-Cdc24 complex by mutating the invariant tryptophan 
(beml'^'^^'^) in the SH3 domain required for binding to proline- 
rich motifs (Larson and Davidson, 2000). beml^'^^^ resulted 
in a slight but significant (P = 0.05 at 50 min) reduction in the 
percentage of polarized cells in the presence of LatA compared 
with the control (Fig. 5, C and D), but biosensor measurements 
showed GEF activity similar to that in the wild-type polar cap 
(P = 0.3; Fig. 4 D). Importantly, the localization of Cdc24 was 
also normal in the beml'^'^''^ background (Fig. 5, E and F). These 
results suggest that GEF localization and activity at the polar cap 
is largely independent of the interaction of Beml with its SH3 
domain ligands. It is interesting to note that beml*'^^'^ mutant 
cells frequently displayed two polar caps in LatA, suggesting that 
the affected interaction is important for the singularity of the 
polar axis in the absence of actin, but further investigation of this 
phenomenon is beyond the scope of this work. 

Symmetry breaking without: actin 
requires Rdi 1 -mediated cytosoiic 
targeting of Cdc4S 

If not through the Beml feedback loop, how then does Cdc42 
polarize without actin? In addition to the membrane-bound 
pool, which is targeted by actin, Cdc42 exists in the cytosol in a 
soluble pool as the Rdil-bound complex (Koch et al., 1997; 
Tiedje et al., 2008; Slaughter et al, 2009). Consistent with a 
new study published recently (Freisinger et al., 2013), rdilA 
cells fail to polarize Cdc42 in the presence of LatA after release 
from pheromone Gl arrest (Fig. 6, A and B), suggesting that 
targeting from the Rdil-bound cytosoiic pool is not only re- 
quired for maintaining but also for establishing Cdc42 protein 
polarization when actin is disrupted. The polarization defect in 
rdilA cells in the presence of LatA was not caused by insuffi- 
cient GEF activity: biosensor measurements indicated the GEF 
activity in rdilA cells to be similar to that in the wild type (P = 
0.89; Fig. 4 D), and expression of Cdc24'^™' did not rescue the 
polarization defect of rdilA cells (Fig. 6, A and B). Furthermore, 
deletion of RDU did not affect growth of rsrlA or axl2A raxlA 
cells, suggesting that Rdil is not required for cellular symmetry 
breaking when actin is intact. However, deletion of RDU exacer- 
bated the slow growth and polarization phenotypes in bemlA 
cells (Fig. 6, A and C). Collectively, the aforementioned results 
show that targeting of Cdc42 from the cytosoiic Rdil-bound pool 
is central to Cdc42 polarization in the absence of actin-based 
vesicular trafficking and that this process is distinct from GEF 
activation, which does not have to be spatially confined. 



Analytical model of actin-independent 
asymmetry breaking without localized 
C3EF activation 

As a theoretical exploration of the polarization mechanism 
based on findings of this study, we used a minimalistic approach 
(Goi-yachev and Pokhilko, 2006; Otsuji et al., 2007; Altschuler 
et al., 2008; Mori et al., 2011) to discover the conditions that 
could allow symmetry breaking through Rdil -mediated Cdc42 
targeting from the cytosol but not local activation of the GEF. 
Main components of the ID model include Cdc42 on the mem- 
brane (local concentration uix.t)), a fraction (ai < 1) of which is 
in the active GTP-bound form, Cdc42 in cytosol v(x,t), free 
(cytosoiic) Rdil R/xj), and the cytosoiic Rdil-Cdc42 complex 
Rc(x,t). Considering the simplest case in which the activity G that 
dissociates the Rdil-Cdc42 complex is proportional to squared 
density of active Cdc42, G{ii) = 02(01")", i-S-, G{u) = Air, in 
which A = a-i_ay\ the membrane targeting term for Cdc42 then 
reads Au^R,. (Fig. 7 A). Cdc42 internalization is given by the 
term fSuRf, in which /3 is the extraction rate of membrane-bound 
Cdc42 by free Rdil. Because cytosoiic Cdc42 exists in Rdil- 
and vesicle-bound states (Das et al., 2012; Slaughter et al., 2013), 
we assume R^ = yv, in which 0 < y < 1. Based on recently pub- 
lished data on the mobile pools of Cdc42 by using fluorescence 
correlation spectroscopy (Das et al., 2012), y was estimated to be 
0.4. Using the conservation condition for mean concentration of 
Rdil (R) in the cell, we have Rf=R-R^ = R-yv. 

The cell is represented as a line segment with a size L. 
Dynamics of the Cdc42 concentrations are described in the re- 
gion 0 < X < L: du/dt = yA(u — utifv — /3(m — ui,){R — yv) + 
D^,{d'uldx^), and dvldt = —yA{u — ui,fv + ji{u — ui,){R — yv) + 
Dy{d\/dx^), in which ut, denotes a nonzero but small basal uni- 
form concentration of Cdc42 in the membrane (see Supplemental 
material). The diffusion of Cdc42 in the cytosol is much faster 
than that in the membrane: » D„. The equations are subject 
to no-flux boundary conditions, and the total amount CL of 
Cdc42 in the cell is conserved. A linear stability analysis showed 
that the growth rate of small perturbations to an initial uniform 
distribution is proportional to the activation level of Cdc42 (see 
Supplemental material). Simulations of the model showed the 
perturbations lead to growth of a single peak of Cdc42 to a steady 
level (Fig. 7 B). 

We explored the parameter space for R (cellular Rdil 
level) and A (essentially a parameter describing the activation 
level of Cdc42 that also impacts Cdc42 targeting) required for 
symmetry breaking. For a fixed value of A, simulations showed 
that formation of stable polarity responds nonmonotonically 
to R, such that polarization occurs when R is in the range 
of 0.5-0.8 (note R values are normalized to the global concen- 
tration of Cdc42, C = 1, and the upper but not lower bound- 
ary depends on the value of A; see Discussion). When R is 
above this range, no polarization occurs as a result of a lack of 
Cdc42 targeting to membrane, whereas when R is below the 
range, a high level of Cdc42 uniformly distributes on the 
membrane (Fig. 7 C). To perform a qualitative experimental 
test of this prediction, we induced expression of Rdil under 
the GALl promoter for varying time periods before release 
from pheromone arrest into LatA-containing media (Fig. S3 A). 
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Figure 5. The Beml-Cla4 interaction is not required for Cdc24 localization or symmetry breaking. (A) Polarization of mCherry-Cdc42 (left] and Cdc24''''°'- 
GFP-Cla4 (rigfit) in RSR 1 BEM 1 cells and in rsrlA bem lA cells in DMSO or LatA upon release from Gl pheromone arrest. Experimental procedure was as 
described in Fig. 1 B. The plots show means from averaging three experiments, and error bars correspond to SEMs. More than 80 cells were counted per 
time point per experiment. (B) Maximum projections of representative cells from A. (C) Polarization of GFP-Cdc42 in axl2A rax TA BEM 1 cells, axl2A rax lA 
berr)TA, and ax/2/1 raxlA bem /"'"^'^ cells, as in Fig. 1 B. (D) Maximum projections of representative polarized cells from C. (E) Maximum projections 
of representative cells with polarized Cdc24-GFP in axl2A rax lA beml^'''^'^ cells after pheromone arrest release in LatA or DMSO. (F) Quantification of 
Cdc24-GFP polarization in axl2A rax lA cells and in axl2A rax lA bem /^"^'^ cells, as in Fig. 1 K. For the comparison of peak values at 0°, P = 0.5. Plots 
show normalized profiles averaged over n > 17 cells. Error bars correspond to SEMs. Bars, 5 pm. 



Quantifying the percentage of polarized cells at 50 min after 
the release, it was apparent that cell polarization occurs opti- 
mally at the concentration of Rdil induced for 60 min with 



galactose and declines sharply above and more gradually below 
this expression level (Fig. 7 D). Measurement of mean fluores- 
cence intensities of Rdil-GFP and GFP-Cdc42, each expressed 
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Figure 6. Rdi1 -mediated cytosolic targeting of Cdc42 is essential for symmetry breaking without actin. (A) Polarization of GFP-Cdc42 in wild-type cells, 
rdilA cells, rdilA cells expressing CDC24^''^^ from the CDC24 promoter, and rdilA bemlA cells in DMSO or LatA after release from Gl pheromone 
arrest. Experimental procedure was as described in Fig. 1 B. The plots show means from averaging three experiments, and error bars correspond to SEMs. 
More than 80 cells were counted per time point per experiment. (B) Maximum projections of representative cells from A. Bar, 5 pm. (C) Serial dilution of 
cells with the indicated genotypes from an overnight culture with an OD of 1 were plated on YPD media and grown at 23°C for 3 d. 



under its endogenous promoter, confirmed that Rdil concen- 
tration is ^0.6-fold of that of Cdc42, within the allowable range 
for symmetry breaking (Fig. S3 B). We also fixed the value 
of R and varied A and found that symmetry breaking requires 
a threshold level of Cdc42 activation (Fig. 7 E), which is qual- 
itatively consistent with the experimental findings in Fig. 3 
(see Discussion). 

□iscussion 

□istiinct: pathways of Cdc4B act:ivat:ian 
and localization contribute to 
symmetry breaking 

The aforementioned results indicate that mechanisms of Cdc42 
localization and activation are not obligatorily coupled. In any 
single location on the cortex, induction of downstream events 
or feedback circuits that rely on Cdc42 effectors depends on 
the concentration of Cdc42'^^^, a product of two factors: (1) total 
Cdc42 protein concentration (controlled by localization on the 



cortex) and (2) relative abundance of the GTP-bound state 
of Cdc42 (controlled by GEF/GAP balance) in that location. 
Modulation of either factor can be envisioned to affect the local 
Cdc42'^^'' level and therefore the capacity for initiating the 
downstream reactions required for symmetry breaking (Fig. 8). 
This concept is sufficient to interpret all the experimental obser- 
vations made in this study. Specifically, our results show that when 
both the actin-based vesicular trafficking and Rdil -dependent 
cytosolic targeting mechanisms are intact, symmetry breaking, 
although less efficient, occurs even without strong GEF local- 
ization or the benefit of Beml's GEF boosting ability (Fig. 8 B). 
Alternatively, optimal GEF activation may compensate for sub- 
optimal protein localization such as in the case of actin disrup- 
tion with LatA (Fig. 8 C). However, when both activation and 
protein targeting are inhibited to a certain degree, such as the 
combination of bemlA and actin inhibition, the threshold of 
Cdc42''^'' concentration could not be achieved to enact the 
downstream reactions and positive feedback loops for symme- 
try breaking (Fig. 8 D). 
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Figure 7. Mafhemafical model of Rdi 1 -dependent symmetry breaking and parameter validation. (A) Schematic of analytical model for Rdil -dependent 
polarization. Distribution of Cdc42 across the cortex is given by u(x,/] (blue curve), of which a constant fraction Oj is activated. Cdc42 extraction is propor- 
tional to free Rdi 1 , Rf (purple crescents), and u[x,fj by a constant /3. Targeting of Cdc42 to the cortex is assumed to be proportional to the squared distribu- 
tion of active Cdc42°^'', (dji;)', the profile of which is represented by the filled orange curve. Cytosolic diffusion D, and cortical diffusion D„ are shown in 
black. (B) Simulation of cell polarization via the Rdil -dependent mechanism. Initial distribution of cortical Cdc42 at the time of perturbation is shown in 
red, with steady-state shown in black and intermediate time points given in blue. The parameter values used in the simulation shown are C = 1 , y = 0.4, 
jS = 0.2, D„ = 0. 1 , D„ = 0.01 , ut = 0.05, L = 3, A = 1 .75, and R = 0.6. (C) Simulated data for Rdi 1 -dependent polarization, showing dependence of 
polarization strength u^^ix on R, cmd Rdil concentration relative to Cdc42. Gray dashed lines indicate cortical Cdc42 concentration in unpolarized states, 
whereas the red curve indicates mean u^^ix over 25 simulations for parameter values in which polarization is allowed. The parameter values are as in B, 
with changing value of R. (D) Experimental assessment of the impact of Rdi 1 expression level on polarization without actin. Expression of Rdi 1 was induced 
under the GAL 1 promoter for the indicated amounts of time (x axis) concurrent with a 1-h Gl pheromone arrest followed by a 0.5-h pheromone arrest in 
glucose media before release into LatA-containing media (see diagram in Fig. S3). The percentages of polarized cells were counted at 50 min after release. 
P = 0.06, comparing polarization in wild type with that of 60-min galactose (Gal) induction. The plots show means from three experiments, and error bars 
correspond to SEMs (blue shading indicates error bars for controls). More than 80 cells were counted per time point per experiment. (E) Simulated data 
for Rdi 1 -dependent polarization, showing dependence of polarization strength u^^^ on A. Gray dashed lines indicate cortical Cdc42 concentration in 
unpolarized states, whereas the red curve indicates mean u^„^ over 25 simulations for parameter values in which polarization is allowed. The parameter 
values are as in B, with R = 0.6 and changing value of A. 
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Figure 8. A schematic explanation of the phenotypic observations based on the cooperation between Cdc42 targeting and activation. The concentra- 
tion of active Cdc42°^'' (filled green circles) on the cortex is a fraction of total Cdc42 (filled and open green circles) determined by the GEF activity of 
Cdc24 boosted by Beml . Cdc42 localization is controlled by both octin-dependent vesicle (blue circles) trafficking and Rdil -dependent (purple crescents) 
pathways. Local concentration of active Cdc42°^'' must reach a threshold level (dotted lines) before imposing sufficient feedback for symmetry breaking. 
(A) V^ild-type cells— both Cdc42 activity and localization ore maximized. The level of Cdc42°^'' is well above the threshold. (B) bem lA cells— GEF activity 
is reduced but with the localization of total Cdc42 remaining high, the threshold level of Cdc42°^'' is still attained. (C) LotA-treoted wild-type 
cells-elimination of the octin-based target pathway reduces localization of Cdc42, but with high activation, the Cdc42°^'' threshold con still be reached. 
(D) LotA-treoted bem lA cells— decreasing both the activation level and the localization of total Cdc42 prevents attainment of the threshold level of Cdc42°^'' 
for symmetry breaking. 



The function of Bem1 in cellular 
symmetry breaking 

The previous conclusion that Beml is essential for symmetry 
breaking was built on the synthetic lethal interaction between 
bemlA and rsrlA and the implicit assumption that rsrlA affects 
only bud site selection but not the core symmetry-breaking 
mechanism. The results presented here and in previous studies 
(Ozbudak et al., 2005; Howell et al., 2012) have shown that 
Rsrl's function is beyond bud site selection but that beml A or 
bemlA rsrlA double mutant cells are capable of polarization 
as long as actin is intact. Our approach of interaction-specific 
perturbation by using point mutations has revealed different, as 
opposed to concerted, functions for the various Beml -mediated 
interactions in cell polarization. The mutant analysis shows that 
the Beml-Cdc42 interaction is important for localization of 



Beml to the polar cap, but this is neither crucial for Cdc24 local- 
ization nor required for symmetry breaking with or without actin. 
This result also implies that Beml can fulfill its main function 
in symmetry breaking without itself being localized. In contrast, 
disruption of the Beml-Cdc24 interaction prevented symmetry 
breaking only when actin was inhibited. This defect correlates 
with significantly reduced GEF activity and can be rescued by a 
Cdc24 construct that enhanced GEF activity but was unable to 
localize. These results strongly suggest that a key function of 
Beml is to boost the GEF activity of Cdc24 globally or locally. 

Beml binding was also thought to enhance GEF activity 
by bridging a complex between Cdc24 with the p21 -activated 
kinase Cla4, thus facilitating Cdc24 hyperphosphorylation. 
Although mutagenesis of as many as 38 phosphorylated sites on 
Cdc24 did not result in any observable phenotype (Wai et al.. 
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2009), disruption of Beml SH3 domain binding to polyproline 
motifs via the beml^'^^'^ mutation resulted in a partial decrease 
in polarization efficiency in the presence of LatA (Fig. 5 C). 
However, as the beml^'^^^ mutation would also disrupt the in- 
teraction with Boil and Boi2, which are together required for 
viability (Bender et al., 1996), the effect of this mutation does not 
permit simple interpretation. Nevertheless, we speculate that, in 
addition to GEF activation, Beml functions as a bona fide adap- 
tor with multiple weak ligand interactions to enhance the affinity 
of protein complexes required for robust cell polarization. 

The emerging role of Rdi1 and Rsr1 
in cellular synnme'fcry breaking 

Our results show that Rdil plays a critical role in actin-independent 
polarization of Cdc42, consistent with a recent study (Freisinger 
et al., 2013). As a Rho GDP dissociation inhibitor, Rdil is the 
cytosolic chaperone of prenylated Cdc42 and thus governs the 
targeting of this pool of Cdc42 to the site of polarization. Our 
previous work showed that Rdil is required for rapid recycling 
of Cdc42 that helps maintain a dynamic polarized Cdc42 con- 
centration in the presence of significant membrane diffusion 
(Slaughter et al., 2009; Das et al., 2012). To recycle Cdc42 back 
to the polar cap, however, the Rdil-Cdc42 complex must be bro- 
ken apart, and what catalyzes this reaction remains a key missing 
link in the cytosolic Cdc42-targeting pathway. We envision that 
to enable symmetry breaking via this cytosolic targeting path- 
way, active Cdc42 on the plasma membrane controls the disso- 
ciation of the Rdil-Cdc42 complex in some manner, restricting 
Cdc42 anchoring to sites of prior Cdc42 accumulation. Model- 
ing of experimentally observed Cdc42 dynamics during steady- 
state polarity predicted that the window of targeting Cdc42 from 
the cytosolic complex must overlap with that of actin-based de- 
livery (Slaughter et al., 2009), but the biochemical mechanism 
underlying this spatial restriction remains unknown. 

The Ras-like GTPase Rsrl has long been known to be es- 
sential for bud site selection, but our results, as well as the finding 
presented in several previous studies (Park et al., 2002; Kozminski 
et al, 2003; Kang et al., 2010), indicate that its role in polarization 
is more extensive and central to symmetry breaking than previ- 
ously thought. Although Rsrl is not strictiy required for polariza- 
tion, especially when actin is intact, in its absence, the established 
polar cap exhibits drastically reduced spatial and temporal stabil- 
ity, as shown in this and previous studies (Ozbudak et al., 2005; 
Howell et al., 2012). Previous studies attributed this instability 
to GAP-mediated negative feedback, and if so, Rsrl may locally 
regulate this negative feedback to enhance the stability of the polar 
cap. The finding that Beml and Rsrl share a role in Cdc24 local- 
ization provides an alternative explanation for the previously re- 
ported synthetic lethality of beml A and rsrl A: perhaps in some 
yeast strain background, without GEF localization, beml A rsrl A 
cells are unable to achieve a sufficient level of cortical Cdc42'^^^ 
to initiate the feedback that brings about symmetry breaking. 

A new model of cytoskelefeon-independent: 
symmetry breaking 

Models of yeast cell polarity have focused on achieving a stable, 
nonuniform distribution of Cdc42 on the membrane (Onsum and 



Rao, 2009). Existing models are mass-conserved reaction- 
diffusion models (Goryachev and Pokhilko, 2006; Otsuji et al., 
2007; Altschuler et al., 2008; Mori et al., 2011). Analysis of these 
models showed that the polarization is caused by Turing-type 
instability (Rubinstein et al., 2012), the physical reason for which 
lies in the significant difference in the Cdc42 diffusion constants 
in the cytosol and in the membrane, in addition to specific non- 
linear dependence of protein recruitment on the local concentra- 
tion of Cdc42 in the membrane. A major feature of the Turing 
instability is that a stable state of the two-component system 
may become unstable in the presence of diffusion. Cdc42 con- 
centrations in the cytosol and on the membrane are considered 
as "master" variables, whereas concentrations of other proteins 
are consequences of the master variable dynamics. 

One model introduced by Goryachev and Pokhilko (2008), 
which formed the basis for subsequent studies (Howell et al., 
2012; Savage et al., 2012), assumes that polarization occurs via 
a mechanism of autocatalytic activation of Cdc42 through re- 
cruitment of the Beml-Cdc24 complex. Rdil was incorporated 
into the model as a passive aide in Cdc42 recycling, prevented 
through a GEF-dependent mechanism from extracting Cdc42 
at the cap. The authors considered a complex model consisting 
of eight reaction-diffusion equations for dynamics of membrane 
bound and cytoplasmic proteins diffusing in 2D and 3D, respec- 
tively (Goryachev and Pokhilko, 2008). Using several simplify- 
ing assumptions, they reduced the original model to a ID version 
for the master variables only. Both the full and the reduced models 
produced qualitatively the same results showing the existence 
of robust clustering of Cdc42 on the membrane. In our approach, 
we focused on the development of the minimalistic ID model 
designed to grasp major features of the actin-independent polar- 
ization process. 

Our model differs mechanistically from that of Goryachev 
and Pokhilko (2008) in the lack of reliance on autocatalytic acti- 
vation of Cdc42 via a proposed Beml-Cdc42-Cdc24 complex. 
Rather, symmetry breaking is achieved through autocatalytic 
Cdc42 protein recruitment from the cytosolic Rdil-Cdc42 com- 
plex. We demonstrate through both model simulation and ex- 
perimental measurements that symmetry breaking depends on 
a window of Rdil concentration relative to that of Cdc42. In 
contrast, cells have more permissive requirements on the level of 
GEF activity, such that although a threshold level is required, 
higher levels are well tolerated. This explains the main experi- 
mental findings that cell polarization can occur with varying 
degrees of GEF concentration in the polar cap and that beml mu- 
tations that impair GEF activation are more detrimental to polar- 
ization when cells are reliant on the Rdil -based targeting pathway 
(i.e., when the actin-based transport pathway is disabled). 

In summary, a key insight resulting from the analyses 
performed in this work is that activation and localization of 
Cdc42 are achieved via distinct mechanisms that contribute 
quantitatively and productively to symmetry breaking. Although 
both aspects of Cdc42 regulation are required, partial deficiency 
in either may be compensated as a result of the presence of 
multiple mechanisms achieving the other This cooperation un- 
derscores both the complexity and robustness of the yeast cell 
polarity system. 
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Materials and methods 

Detailed model description 

We considered a 1 D model describing Cdc42 protein dynamics in a yeast 
cell undergoing symmetry breaking transition from a uniform state to a po- 
larized state. For the actin-independent pathway, we introduce the follow- 
ing components of the model: Cdc42 on the membrane (local concentration 
u(x,/)] with active fraction oi < 1 , Cdc42 in cytosol v(x,t), free (cytosolic) 
Rdil Ri{x,fj, protein complex Rdil-Cdc42 Rc(x,f), and a certain activity 
G[u] leading to breaking up of the complex into free Rdil and membrane- 
anchored Cdc42. We considered the simplest case when this activity is 
proportional to squared density of active membrane Cdc42, i.e., G(u] = 
a2ai^u^ = Au^, in which A = 020|^. 

The membrane targeting term reads G{u]Rc, and the internalization 
term is fiuRf, in which /3 is the extraction rate. Assuming that the Rdil- 
Cdc42 complex exists as a fraction of total cytosolic Cdc42, we find Rc = 
■yv, in which 0 < y <1 . Given conservation of total amount of Rdil [R] in 
the cell, we also hove Rf = R — Rc = R — yv. This means that both Rc{x,t] and 
Rf(x,f) are dependent on the Cdc42 dynamics. 

Spatiotemporal dynamics of Cdc42 is described by the equations in 
the region 0 < x< L: 



du 
dt ' 



yAu^v - 



Pu{R - yv] + 



'5x2' 



5^ x 2 
— = -yAu V 



- Pu[R - yv] + D^ 



d^v 
'5x2' 



(1) 



in which the diffusion of the Cdc42 cytosolic form is much faster than the 
membrane one: D, >> D„. The cubic nonlineority is shown to be critical for 
the Turing-type instability in polarization models (Goryachev and Pokhiiko, 
2008). The system of Eq. 1 is subject to no-flux boundary conditions. Sum- 
ming up the equations and integrating over the spatial interval, we obtain 
the Cdc42 mass conservation condition 



|q(u + v)c/x = CL = constant. 



(2) 



in which the constant C is the model parameter representing the mean 
Cdc42 concentration. The reaction terms used in Eq. 1 lead to a basic uni- 
form solution u = 0 corresponding to absence of Cdc42 on the membrane. 
To have some nonzero basal level ut on the membrane, we modify the re- 
action term in Eq. 1 to arrive at 



du 
df ' 



dv_ 
dt ' 



f(u,v)-^D„ 



5x 
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5x^ 
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(3) 



Being stationary and spatially uniform, the solution {uo,Vo} is independent 
of time and spatial variables, so that all derivatives vanish, and this solu- 
tion should satisfy the equation f(uo,Vo) = 0. From Eq. 2, it follows that the 
basic solution satisfies the condition uo -i- vo = C, and we find from Eq. 3 



"0 



^0 



AC - p + Au^ - s 
2A ' 

AC + P - Au{, -I- s 
2A ' 



..^I^AC + p-Aui,f -ApAR/y. 



(4) 



As the uniform steady-state concentrations should be nonnegative, we find 
a necessary condition on the model parameters R > yC, which was shown 
to be satisfied with experimental measurements and the estimated value of 
0.4 for y. 



A previous study (Rubinstein et al., 201 2) showed that the maximal 
linear growth rate a-„ of the perturbation can be approximated by the par- 
tial derivative of f with respect to the variable u: cr^ ~ df/du = 2yA(uo — 
^b]vo ~ I3{R — yvo), computed at the basic solution (Eq. 4). Meanwhile, the 
basic solution satisfies the relation yA(uo — Ubjvo = y8(/? — yvo), which pro- 
duces on estimate for the growth rate for membrane Cdc42 in the form 
o-„ = p{R - yvo] = yA{uo - l/Jvo. 

The last relation shows that the growth rate of membrane Cdc42, 
which determines the kinetics of polarization establishment, is proportional 
to the concentration of Rdil [R). It also depends on the activation level A, 
and more precisely, it is determined by the balance between the activation 
level A and the Rdil expression R. 

For a fixed value of activation, with increase in R, the internalization 
term wins over the membrane-targeting term, eventually blocking polarization. 
On the other hand, for small R < yC, polarization is also not possible. This 
means that there exists a range of the Rdi 1 expression level R for which the 
polarization is possible, and the upper boundary of this range varies with A. 
This behavior is illustrated in Fig. 7 C in which the maximal value u„ox of the 
membrane-bound Cdc42 concentration obtained by simulation of Eq. 3 is 
shown for a fixed value of A and increasing R. The nonpolarized state (with 
the ratio u^ax/ Um/n < 1 .2) is shown by the dashed curve. The parameter val- 
ues used in this simulation are C = 1 , y = 0.4, /3 = 0.2, = 0. 1 , D„ = 0.01 , 
ut = 0.05, and 1=3. The simulations for each set of the parameters were 
performed 25 times, and the mean values of u^„^ were computed. 

Genetic manipulations 

Site-directed mutagenesis was performed using the site-directed mutagene- 
sis kit (QuikChange II XL; Agilent Technologies), and the final product was 
sequenced to ensure that there were no secondary mutations introduced. 
All yeast strains used in this study ore described in Table SI . Techniques 
for yeast cell culture and genetics were essentially as previously described 
(Burke et oL, 2000). Transformation of plasmid DNA into yeast was per- 
formed based on the lithium acetate method (Ito et oL, 1 983). Transforma- 
tion of DNA fragments was performed by the some method, but after 
transformation, cells were recovered in nonselective medio for at least two 
cell cycle times before plating on selective medio, and genomic integration 
was confirmed by PGR. The plasmid DLB3170 was a gift from D. Lew 
(Duke University School of Medicine, Durham, NC). 

a-Factor release assays 

Cells were arrested at 23°C for 1.5 h using 20 [Jg/ml a-foctor and re- 
leased into the cell cycle by washing three times in sterile water before re- 
suspending in fresh media with either 50 pM LatA or equivalent DMSO. 
LotA treatment was confirmed to result in actin polymerization within 
15 min (Fig. SI B). Samples containing Beml-GFP or Cdc24-GFP were 
imaged as live cells 50 min after release. Samples were token at 20, 50, 
and 80 min after release, fixed in 4% paraformaldehyde for 1 5 min before 
washing in PBS, and stored at 4°C for <48 h before imaging. For each 
time point, >80 cells were scored for polarity. The assay was repeated at 
least three times for each genotype. For release assays with Rdil expres- 
sion induced under the Gall promoter, pGALT-RDIl cells were grown 
overnight in 4% raffinose media, to which 4% galactose was added at the 
appropriate time point relative to the addition of a-foctor such that addition 
of glucose 30 min before release from a-foctor would end the induction. 

Microscopy and live-cell imaging 

Imaging was performed at 23°C on a spinning-disk confocal microscope 
(UltroVIEW; PerkinElmer), including on inverted microscope (Axiovert 
200 M; Carl Zeiss), attached to a spinning-disk confocal system (CSU-Xl ; 
Yokogawo Corporation of America) and electron multiplying charge-coupled 
device (EM-CCD) camera (C9100; Homomotsu Photonics) with Volocity 
acquisition software (PerkinElmer] or a similar system with MetaMorph ac- 
quisition software (Molecular Devices). Single time point images of live or 
fixed cells were collected as a series of optical sections, with a step size of 
0.5 pm. ImageJ software (v. 1 .37; National Institutes of Health) was used 
to process the images. Final images are maximum projections that have 
been background subtracted and contrast adjusted for clarity. Cdc24-GFP, 
Beml-GFP, and mCherry-Cdc42 images were token with a lOOx aPlon 
Fluor, NA 1 .46 objective lens, whereas GFP-Cdc42 and Cdc24-GFP-Cla4 
images were token with a 63x Plan Apochromot, NA 1 .4 objective lens. 

For live-cell videos of GFP-Cdc42, cells were arrested in Gl via 
pheromone for 1 .5 h and released into LotA-contoining synthetic complete 
(SC) media for 30 min before slide preparation on a 1% agarose and 
1 00 pM LatA pad. Z-stock images were acquired at 30-s time intervals with 
0.7 pm between slices. Maximum z projections of each time point were 
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generated, cells showing polar cap formation and disappearance within a 
single movie were identified, and the maximum intensity of the polar cop 
was measured and plotted over time. These were smoothed via a rolling 
mean of six time points and then fitted to Gaussian curves in OriginPro 
(OriginLab]. The duration of the cap is reported as the full width half maxi- 
mum of the Gaussian fit. 

For live videos of Cdc24-GFP, cycling mid-log phase cells were im- 
mobilized on ]% agarose in SC media and imaged at 2-min time intervals 
as z-stack images with 0.7 pm between slices. Cells were aligned before 
kymograph generation by manually tracking either the polar cap or bud 
neck after bud emergence. Kymographs were then generated by averaging 
over a thickness of 3 pixels along the cortex using a custom plugin in ImageJ. 
This and other plugins can be found at the Stowers Institute ImageJ Plugins 
website and in a ZIP file provided online in Supplemental material. 

Bern 1 polarity analysis 

Maximum projections of Beml-GFP fluorescence images from the experi- 
ment in Fig. 2 I were background subtracted and processed using the 
ImageJ Smooth function to remove noise. Unbudded cells were circled, 
and the maximum and mean pixel intensities were calculated for a total of 
520 cells per genotype. The Bem 1 polarity was calculated as the maximum 
pixel intensity divided by the mean intensity per each cell. 

Cdc24-GFP and mCherry-Cdc42 profile analysis 

Two-channel z-stack images of cells expressing Cdc24-GFP and mCherry- 
Cdc42 were background subtracted, and maximum projections were gen- 
erated. A cortex mask of the cell was then generated from the mCherry-Cdc42 
channel using a custom macro based on the Li Dark thresholding algorithm 
in ImageJ (Li and Tam, 1998). The mask was applied to the Cdc24-GFP 
channel, and the mean intensity of the masked region in 10° increments 
around the cortex was calculated using a custom plug-in to generate a corti- 
cal intensity profile or Cdc24-GFP or similarly with the mCherry channel. 
This was completed for >15 cells, including all visibly polarized cells from 
each of several images. The profiles were aligned by maximum value, and 
each was normalized to the mean intensity of the cortex in the 1 80° region 
opposite the peak. The normalized curves were averaged. Statistical analy- 
sis was performed on peak (0°) values. 

Cdc42 activation biosensor imaging and analysis 

Cells were grown to mid-log phase before analysis. Polarized cells were im- 
aged on a spinning-disk confocal microscope (UltraVIEW) including an in- 
verted microscope (Axiovert 200 M), a spinning-disc confocal system 
(CSU-Xl), a charge-coupled device (2x2 bin; ORCA-R2; Homamatsu Pho- 
tonics), and a PhotoKinesis accessory (PerkinElmer) using Volocity acquisition 
software with a 1 OOx aPlan Fluar, NA 1 .46 objective. Cells were imaged at 
maximum speed for a total of 42 frames, with 488-nm laser power at 7% 
and exposure at 100-300 ms depending on biosensor expression level per 
cell. mCherry throughout the cell was bleached after frame 20, using one 
iteration of 568-nm laser at 1 00%. mCherry fluorescence was checked after 
acquisition to ensure complete bleaching. Controls indicated no significant 
bleaching of GFP by 568-nm excitation nor cross talk of mCherry fluores- 
cence into the GFP channel. The sensor was expressed within a similar range 
(measured by fluorescence) in all strains, and cortical fluorescence intensity 
did not correlate with FRET (Fig. S2). Analysis was performed using ImageJ 
software. The first two frames after bleaching (frames 21-22) were discarded 
because of a potential delay of 568-nm laser shut off. The sum of the full time 
series was then used to create a cortex mask of the cell, using a custom 
macro based on the Li Dark thresholding algorithm in ImageJ (see Supple- 
mental material; Li and Tarn, 1 998). The mask was applied to the series, and 
the mean intensity of the masked region in 1 0° increments around the cortex 
was calculated at each time point using a custom plug-in. The mean intensity 
and the SD intensity of each 1 0° increment before and after bleach were cal- 
culated in RStudio (RStudio, Inc.). The SDs of all 10° increments before and 
after bleach were averaged, and cells whose mean SD was >10% of its 
mean cortical fluorescence were discarded. Intensity profiles of remaining 
cells were then aligned such that polar cap center was located at 0°. FRET ef- 
ficiency was calculated for each 10° segment in remaining cells as 100 x 
(postbleach - prebleach)/prebleach. FRET efficiencies for each 10° incre- 
ment were then averaged over >25 cells, as in the plot in Fig. 4 C. The peak 
FRET efficiency (Fig. 4 D) was determined as the mean of the 30° region con- 
taining the 10° increment of highest mean FRET efficiency among aligned 
profiles, as illustrated by the yellow box in Fig. 4 C. 

Statistical analysis 

Statistical differences between two sets of data were analyzed with a two- 
tailed unpaired Student's ttest. 



Online supplemental material 

Fig. SI shows additional data for characterization of spatial cue indepen- 
dent polarization backgrounds, rsrTA and axl2A raxlA. Fig. S2 shows 
mean cortical intensity of the Cdc42 biosensor plotted against FRET effi- 
ciency for each cell in the analysis. Fig. S3 shows a schematic of the ex- 
perimental procedure used to generate Fig. 7 D and the relative levels of 
Rdil and Cdc42. Videos 1 and 2 show polar cap dynamics of GFP-Cdc42 
in LatA in an rsrlA cell and on axl2A raxlA cell, respectively. Videos 3 
and 4 show localization of Cdc24-GFP during budding in rsrTA bemlA 
cells and axl2A raxlA bemlA cells, respectively. Table SI shows yeast 
strains used in this study. A ZIP file is also provided that contains custom 
plug-ins and macros written for ImageJ, used to calculate cortical intensity 
profiles and kymographs. Online supplemental material is available at 
http://www.jcb.org/cgi/content/full/jcb.201 3041 80/DCl . 
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